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Current Thrusts 


• BACKGROUND 

Decrease design cycle time => Rapid turn-around 

Increase design/process fidelity => High accuracy and low variation 

(Accurate flow physics) 

Increase discipline integration => Increased range of options via IT 

• OBJECTIVE 

Develop advanced computational methods to provide appropriate fidelity 
computational analysis/design capability 



Arne* Re»mh Cm'tt 


Current Thrusts 


• APPROACH 

• Advanced computational methods 

Enhance / accelerate viscous flow simulation procedures 
Develop hybrid/polyhedral-grid procedure for viscous flow 

- Methods for evaluation and design process 

Develop real time transonic flow simulation procedure for production 
wind tunnel 

Develop intelligent data management technology (IT) 

- Validation of methods and flow physics study 
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Recent Accomplishments / Work In Progress 

• VALIDATION OF METHODS AND FLOW PHYSICS STUDY 

- Wingtip and wake vortex study 

- Unsteady incompressible N-S solver and adaptive control 

- Aeroacoustic applications : 

Flap edge noise and acoustic modeling 
DNS/LES for jet aeroacoustics 
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Computational Methods Development 


TECHNICAL HIGHLIGHTS 


Acceleration of Viscous Flow Simulation 


• FUTURE VISION / PAYOFFS / BARRIERS 

- Two orders of magnitude improvement in turn around time will enable RANS 
applications to viscous flow design problems 

- Improved proc ess time will impact design cycle time 

- Existing Navier-Stokes solution procedures are not efficient enough for 
engineering applications. 

• APPROACH 

- Tum-around time of current RANS solver can be reduced by advanced 
convergence acceleration algorithms combined with parallel implementation 

- Advanced methods will be developed to reduce CFD process time from 
CAD to solution 

• TASKS 

- Enhancement o ,f OVERFLOW : preconditioning & multigrid, MPI version 

- Development of Cartesian-based hybrid-grid method 

- Development of acceleration algorithms 



Acceleration of Viscous Flow Simulations 


Convergence Acceleration for N-S Computations 
NACA4412 M = 0.001 o=fl Re=1.5 Million 
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• Convergence Acceleration of OVERFLOW 

DESCRIPTION 

Convergence to steady state of the compressible 
Navier-Stokes code, OVERFLOW, has been 
accelerated with multigrid and also with low-Mach 
number preconditioning. 

ACCOMPLISHMENTS 

A multigrid algorithm has been added to OVERFLOW. 
It handles “Chimera" geometries. 

A low Mach number preconditioning option has been 
added for cases with low free stream Mach numbers 
(less than 0.2). 

SIGNIFICANCE 

Convergence to steady state has been accelerated by 
a factor of upto 6 for NACA 4412 at M=0.2 and Re=1 .5 
million case. 

WORK IN PROGRESS / FUTURE 

Robustness of the algorithms in multi-zone geometries 
is being enhanced. 

POC:D. Jespersen, 4-6742, jespersenOnas.nasa.gov 
T. Pulliam, 4-6417, pulliamOnas.nasa.gov 






Acceleration of Viscous Flow Simulations 


Effect of Low Mach Number Preconditioning in 
OVERFLOW 

NACA4412 a=0.0 Re=1 5x10® Spala rt- Alima ras 
NACA4412 o-O.O flt-li Miion Spalut-AJmaru 
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• OVERFLOW /MPI 

DESCRIPTION 

To reduce tum-around time for obtaining Navier- 
Stokes Solutions, a parallel version of OVERFLOW 
has been developed that exploits parallelism at all 
levels. 

ACCOMPLISHMENTS 

An MPI version of OVERFLOW has been 
successfully run on workstations, SP2, T3E, Power 
Challenge. 0rigin2000 and On yx2. 

SIGNIFICANCE 

Distributed linear solver is at 60-80% of ideal 
efficiency. 

Load balancing allows splitting of large zones across 
processors and clustering of multiple small zones 
onto one processor. 

WORK IN PROGRESS / FUTURE 

The code will be integrated into the standard 
OVERFLOW and will be a part of the OVERFLOW 
distribution. 

POC: D. Jespersen, 4-6742, jespersenOnas.nasa.gov 
T. Pulliam, 4-6417, pulliamOnas.nasa.gov 



Acceleration of Viscous Flow Simulations 


Illustration of Polyhedral Interface Cells 



o Finite Volume Methods for Arbitrary Polyhedral 
Grids 

DESCRIPTION 

One key element in developing a hybrid-grid method is 
how to handle irregular polyhedral cells at interface 
boundaries. The objective of the current effort is to 
develop theoretical foundation to achieve very high 
order accurate solutions for general polyhedral grids. 
ACCOMPLISHMENTS 

High order flow variables and fluxes for arbitrary 
polyhedral and polygonal grids have been formulated 
based on exact integration of polynomials and 
quadrature approximations. 

SIGNIFICANCE 

Processing time to generate grids for flows past 
complex geometries can be significantly reduced by 
relaxing grid generation requirements, i.e. structured 
viscous grids on Cartesian outer grids. 

FUTURE PLAN 

Formulas will be incorporated into existing RANS code 
to generate polyhedral grid Navier Stokes capability. 
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POC: Y. Liu / M. Vinokur, 4-6667, iiuGnas.nasa.gov 



Acceleration of Viscous Flow Simulations 


Structured-Cartesian Hybrid Grid Solution 
NLR 7301, M=0.2, Re=10 a , oc=0 



• Hybrid-Grid Method for CAD to RANS Solutions 

DESCRIPTION 

For 3-D applications: Multiblock grid generation 
procedure is expensive; Overset grid approach suffers 
from interpolation inaccuracies; Existing unstructured grid 
method using triangles needs improvement for viscous 
flows. 

The objective is to develop a hybrid grid method using 
locally nonoverlapping structured grids on the Cartesian 
grid background. This requires development of advanced 
methods for the interface boundaries filled with irregular 
polygons. 

ACCOMPLISHED 

High order fluxes for irregular polygons have been 
formulated. Quadratic reconstruction formula has been 
implemented in 2-D with Roe splitting and an implicit 
GMRES method. 

SIGNIFICANCE 

The present approach will simplify the grid generation 
process substantially. 

WORK IN PROGRESS / FUTURE 

High order method for general polyhedral cells will be 
implemented in 3-D. 

POC:M.Delanaye, 4-5176, delanaye© nas.nasa.gov 


Real Time Transonic Flow Simulation in Wind Tunnel 


• FUTURE VISION / PAYOFFS / BARRIERS 

- Cost / solution time determines the hierarchy of method selected 

- Quick solutions can be generated in real time to provide critical information 
for early design/lest decision 

- RANS approach is impractical for real time applications during production wind tunnel 
tests 

• APPROACH 

- For real time applications, implement fast potential solver (TOPS) using the same grid 
topology used for OVERFLOW 

- Develop/implement boundary layer correction 

- Develop zonal methods through IT based decision process 

• TASKS 

- Transonic overset-based potential solver (TOPS) 

Develop TOPS to include wind tunnel wall and support 
Add boundary layer correction into TOPS 
Determine range of applicability and variations 
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Real Time Flow Simulations 


Mach contours: three-zone RAE Wing A with B2 
fuselage at M=0.9, a=<) 


• Transonic Overset Potential Solver (TOPS) 

DESCRIPTION 



Potential/Boundary-Layer approach for transonic cruise 
analysis and design can be up to 100 times faster than 
Navier-Stokes. Full potential capability using existing 
Chimera zonal grid approach will be valuable for real time 
application of CFD for wind tunnel tests. 
ACCOMPLISHMENTS 

Transonic 3-D two-zone wing solver has been completed. 

Transonic 3-D three-zone wing/body solver has been 
completed. 

New doner cell search algorithm developed which is upto 
xl 00 faster than existing algorithm. 

SIGNIFICANCE 

Results to date demonstrate viability of full potential 
Chimera approach, which is compatible with existing 
OVERFLOW environment. 

WORK IN PROGRESS / FUTURE 

Complex geometry applications are in progress. The 
method will be extended to viscous flow, W 7T wall 
interference study and design optimization. 



POC: T. Holst, 4-6032, tholstOmail.arc.nasa.gov 


Data Management and Feature Extraction Method 


• FUTURE VISION / PAYOFFS / BARRIERS 

- Intelligent information management is essential to IT based computational 
simulation methods development 

- Method for data compression and reconstruction within specified accuracy can 
reduce storage requirement and transmission time significantly (i.e. at least an 
order of magnitude reduction) 

- CFD data compression tool (code) is not available 

• APPROACH 

- Develop multiresolution algorithm to compress 3-D vector quantities and 
demonstrate the capability using CFD data such as wing tip vortex flow 

• TASKS 

- Develop data compression and feature extraction algorithm / code 


Arnei Renirth Ccntar 


Data Management and Feature Extraction Methods 


Compression of Analytic Data Using 
Multiresolution Algorithm 
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Compression 
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• Multiresolution Algorithm for 3-D Data 
Compression 

DESCRIPTION 

Numerical solution is, in general, only an approximation 
to the physical solution therefore for many purposes 
truncation of numerical data is acceptable. An efficient 
data compression algorithm for large 3-D data sets is 
desired for storage and transmission of data. 
ACCOMPLISHMENTS 

A prototype 3-D truncation compression algorithm has 
been programmed in Matlab and Fortran. This code is 
based on Marten’s multiresolution algorithm with smooth 
basis functions. 

SIGNIFICANCE 

Since most numerical flow solutions are smooth almost 
everywhere, this algorithm should provide efficient CFD 
data compression. 

Truncation even with high accuracy (e.g. 99%) may 
allow significant compression thus reducing storage and 
transmission requirements 

WORK IN PROGRESS / FUTURE 

Will develop a 3-D CFD data compression code 


CFD Application to Safety 


Wake of a Lifting Wing !"*J^ ^ *t j — « 
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• Characterization of Wake Vortex 

DESCRIPTION 

To develop predictive capability for aircraft wake, it is 
important to characterize the initial tip vortex 
formation/roll-up process and the near wake 
behavior. The objective of the present study is to 
assess current capabilities and define requirements 
for numerical simulations. 

ACCOMPLISHED 

Tip vortex formation and roll-up process have been 
studied in detail addressing the roles of numerics and 
turbulence model in prediction. 

Applicability of a simplified Navier-Stokes solver 
using velocity-vorticity formulation was investigated 
for wake modeling. This model was compared with 
0.03 scale model test of B747 (80x120 data by 
Rossow et al.). 

SIGNIFICANCE 

The current work is an attempt to develop a physics- 
based wake vortex prediction model. 

WORK IN PROGRESS / FUTURE 

Grid adaption methods are being investigated to 
reduce the grid density requirement. 

POC: J. Dacles-Mariani. 4-5369.daclesCnas.nasa.gov 




Wake of a Lifting Wing 


Initial Vortex Formation 

• less than 5 wingspans behind wing 

• tip vortex starts its roll-up process from wihg 



• greater than a few hundred wingspans behind wing 

• diffusion across the vortex takes place 


COMPUTED RESULTS 

Test Conditions for B-747 Modeled in the 80 x 120 ft Wind IXinnel*: 

Frees tream Velocity =131 ft/sec 
Dynamic Pressure = 20 lb/ft2 
Distances = 81 and 162 ft 


Data on Wake Generating Model 


Parameter 

B-747 

Scale of Model 

0.03 

b (wingspan) 

70.5 in. 

c (wing chordlength) 

10.1 in 

Re 

660,000 

Aspect Ratio 

6.96 

S (wing plan form area) 

4.94 ft 2 


Rossow, et ai. (1993) 



WINGTIP VORTEX VALIDATION STUDY 

INS3D- i/P CODE. NACA0012. Re = 4.6e6. a = 10°, GRID SIZi: = 115x137*83 
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WEVGTIP VORTEX STUDY 

Turbulence Model/Grid Density Assessment 





Velocky Magnitude 


COMPUTED RESULTS 


FULL NAVIER-STOKES FORMULATION 


• II A (Campwd) 

• 121.3 ft (Compuietf) 

• 142 ft (Compoted) 


Maximum Vfthie (Mawwed) 

-**- Time Averand Vriae ( Mow e d ) 
- o- • Minimum Vahte (Mea w red) 

Computed 



Computed w- velocity progression in 
the stream wise direction 


Comparison of the up- and downwash 
velocity profile with measured data at 
x/c = 162 ft 


AXIAL PROGRESSION OF WINGTIP VORTEX IN NEAR-WAKE 



NACA 0012 WING 
WITH ROUND TIP 

Re = 4.6 Million 

a= 1(T 
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Incompressible Navier-Stokes Methods 


Axial progression of wingtip vortex in near-wake 
NACA 0012 wing with round tip:Re=4.6x10 a , a=I0° 





• Fractional-Step Method for Time-Accurate 
Incompressible Navier-Stokes Computations 

DESCRIPTION 

Both artificial compressibility method and preconditioned 
compressible N-S methods require subrterations for 
time-accurate solutions of incompressible Navier-Stokes 
equations. 

The objective is to develop an efficient unsteady 
incompressible method for 3-D applications using 
fractional step method. 

ACCOMPLISHMENTS 

Implicit Runga-Kutta scheme of order N and GMRES 
have been implemented. 

Spalart-Allmaras as well as Baldwin-Barth model have 
been implemented. 

Validation cases include Pulsatile flow through 
constricted channel and a wake vortex problem. 

SIGNIFICANCE 

The CPU time requirement for unsteady calculations has 
been reduced by a factor of 3-5. 

WORK IN PROGRESS / FUTURE 

Robustness at high Reynolds number with highly 
stretched grid cases is being enhanced. 


Exploratory Research - Noise Correlation 


Far-Fieid Noise from A Family of Single Element 
NACA Wings 


• Correlation of Wingtip Noise with Mean Flow 
Parameters 



» JO 30 40 
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DESCRIPTION 

Correlations between steady, mean-flow parameters and 
radiated farfield noise were found for a single-element 
wing with a round tip. 

APPROACH 

Steady, incompressible RANS (INS3D-UP) calculations 
provided detailed flow information for the wing 
investigated. An experimental test performed at LaRC 
(Brooks and Marcolini, 1986) provided farfield noise 
measurements. The two investigations were combined to 
generate a correlation that held for combinations of 
aspect ratio, Reynolds number, and angle of attack. 
SIGNIFICANCE 

Correlations help understand what flow features produce 
farfield noise, suggest ways to reduce noise, and make 
possible an acoustic model that would use steady flow 
parameters. 

WORK IN PROGRESS / FUTURE 

The same correlations are being applied to the flap side 
edge of a multi-element wing. 
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POC; D. Mathias, 4-0836, mathias ©nas.nasa.gov 




(51 % of blade height) 


Impeller exit radius : 
5.5 in. 


V* : Impeller exit wheel 
speed (303.5 ft/sec) 







Exploratory Research-CFD Application to Biofluid 



• CFD-Assisted Development of a Long 
TermVentricular Assist Device (VAD) 
DESCRIPTION 

Short-term VAD (two-week pump) has been successfully 
developed with the aid of INS3D. For long term use, the 
device has to be controllable adjusting to the patients 
activity which requires full unsteady flow analysis. 

The objective is to develop a long term VAD by 
implementing adaptive control in conjunction with 
unsteady flow simulations. 

ACCOMPLISHMENTS 

Developed unsteady flow simulation procedure which can 
be used for the entire configuration. 

Pulsatile flow through various inlet cannula designs were 
compared. 

SIGNIFICANCE 

Development of the device could not have been possible 
without CFD capability. 

WORK IN PROGRESS / FUTURE 
Unsteady flow simulation for the entire geometry 
including an adaptive control model will be attempted. 
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POC: 


C. Kiris, 4-4485, kirisOnas.nasa.gov 

D. Kwak, 4-6743, (fltwakOmail.arc.nasa.gov 


Historical Examples 
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• MAJOR ACCOMPLISHMENTS 

Ames has produced the original basic advances in CFD necessary for 
aeronautic applications and has helped extend CFD into many new and 
critical aerospace applications. 

Following are only a few examples in critical CFD processes illustrating 
Ames contribution in the recent past: 

Grid Generation advances and codes: 

Developed elliptic, hyperbolic, unstructured point-insertion and local 
optimization, patched and overset (Chimera) domain decomposition 
grid procedures. 

Resulting codes include 

3D GRAPE, HYPGEN, DCF, SAGE, SURGRID, PRISIM, DELAUNAY 
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Historical Examples 


• MAJOR ACCOMPLISHMENTS (Cont'd) 

CFD Algorithm advances: 

Developed following algorithms and disseminated throughout the US: 

MacCormack explicit algorithm 

Beam-Warming approximate factorization algorithm 

AF1/AF2 algorithm for small disturbance and potential flow 

Pulliam-Chaussee diagonalized ADI algorithm 

Steger-Warming flux-splitting algorithm 

Pseudo compressibility algorithm for Incompressible Navier-Stokes equations 
Linear-reconstruction unstructured scheme 


Flow Solvers and codes: 

Small-disturbance, potential and BL: 


Euler/Navier-Stokes: 

Incompressible: 

Aero-elastic: 

Rotorcraft: 

Turbomachinery: 


Bailey-Ballhaus, LTRAN2JAIR/ 
TWING, BL3D 
ARC2D/3D, TNS, PNS, TIGER, OVERFLOW, CNS 
INS3D family of codes 
ENSAERO 

TURNS, O VE R FLOW_Rtrcrft 
ROTOR, STAGE 


& 
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Historical Examples 


• MAJOR ACCOMPLISHMENTS 

Ames has produced the original basic advances in CFD necessary for 
aeronautic applications and has helped extend CFD into many new and 
critical aerospace applications. 


Following are only a few examples in critical CFD processes illustrating 
Ames contribution in the recent past: 

Grid Generation advances and codes: 

Developed elliptic, hyperbolic, unstructured point-insertion and local 
optimization, patched and overset (Chimera) domain decomposition 
grid procedures. 

Resulting codes include 

3D GRAPE, HYPGEN, DCF, SAGE, SURGRID, PRISIM, DELAUNAY 
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Historical Examples 



• MAJOR ACCOMPLISHMENTS (Cont'd) 

CFD Algorithm advances: 


Developed following algorithms and disseminated throughout the US: 

MacCormack explicit algorithm 

Beam-Warming approximate factorization algorithm 

AF1/AF2 algorithm for small disturbance and potential flow 

Pulliam-Chaussee diagonalized ADI algorithm 

Steger-Warming flux-splitting algorithm 

Pseudo compressibility algorithm for Incompressible Navier-Stokes equations 
Linear-reconstruction unstructured scheme 

Flow Solvers and codes: 


Small-disturbance, potential and BL: 


Euler/Navier-Stokes: 

Incompressible: 

Aero-elastic: 

Rotorcraft: 

Turbomachinery: 


Bailey-Ballhaus, LTRAN2.TAIR/ 
TWING, BL3D 
ARC2D/3D, TNS, PNS, TIGER, OVERFLOW, CNS 
INS3D family of codes 
ENSAERO 

TURNS, OVER FLO W_Rtrcrft 
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• MAJOR APPLICATIONS 

Aerospace: 

- Coupled aerodynamic-structured prediction of flutter 

- Full fighter aircraft performance predictions: transonic cruise (F-16), 
high alpha (F-18 with Dryden), STOVL (Harrier) 

- Full rotorcraft performance prediction (e.g. V-22) 

- Full simulation of the Space Shuttle stack performance (with JSC) 

- "Tip-to-tail" hypersonic aircraft performance prediction 

- Single- and multi-stage turbomachinery performance prediction 

- Simulation/redesign of the SSME Hot Gas Manifold 

- Coupled aero-optics-acoustics analysis of a space observation system,SOFIA 
Non-aerospce: EExamples of technology reapplications 

- First simulation of Artificial Heart (Penn State) 

- Development of the NASA/DeBakey Ventricular Assist Device 

- Simulation of naval vehicles (submarine and propeller with the navy) 



Future Directions 


• RESEARCH 


Algorithm : 

Convergence Acceleration 

Viscous Separated / Unsteady Flow Methods 

Solution Adaptive Method 

Parallel Computing and Real-Time CFD 

High-speed and chemically reacting flow 

LES/DNS (Leveraged by CTR’s Advances) 

Flow Physics: 

Transition 

Turbulence 

Combustion/ Real gas effects 
Aeroacoustics 

IT: 

Computational Tool for Non-experts (User-Friendly) 

Software : 

Outsource (Commercial Vendors) 

Amci Research Center 





